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Abstract: We derive a set of design guidelines and a figure of merit 

to aid the engineering process of on-chip waveguides for strong Stim¬ 
ulated Brillouin Scattering (SBS). To this end, we examine the impact 
of several types of loss on the total amplification of the Stokes wave 
that can be achieved via SBS. We account for linear loss and nonlinear 
loss of third order (two-photon absorption, 2PA) and fifth order, most 
notably 2PA-induced free carrier absorption (FCA). From this, we derive 
an upper bound for the output power of continuous-wave Brillouin-lasers 
and show that the optimal operating conditions and maximal realisable 
Stokes amplification of any given waveguide structure are determined by a 
dimensionless parameter & involving the SBS-gain and all loss parameters. 
We provide simple expressions for optimal pump power, waveguide length 
and realisable amplification and demonstrate their utility in two example 
systems. Notably, we find that 2PA-induced FCA is a serious limitation to 
SBS in silicon and germanium for wavelengths shorter than 2200nm and 
3600 nm, respectively. In contrast, three-photon absorption is of no practical 
significance. 
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1. Introduction 

Stimulated Brillouin Scattering (SBS) is the coherent and self-amplifying interaction between 
light waves and a hypersonic acoustic wave that are confined in the same waveguide or bulk 
material [1]. Initially predicted by Brillouin [2] and first observed in quartz [3], it has been well 
known for many years in optical fibers as a highly resonant and strong third-order nonlinearity 
in both backward [4] and forward directions [5]. More recently, there has been an intense re¬ 
search effort to generate SBS in on-chip waveguides: SBS has been observed in chalcogenide 
glass rib waveguides [6] and, inspired by theoretical considerations [7], silicon nanowires [8,9]. 
These experiments have greatly broadened the applicability of SBS to include a large number 
of on-chip applications, such as powerful narrow-band light sources [10,11], non-reciprocal 
light propagation [12,13], slow light [14] and signal processing in the context of microwave 
photonics [15-17]. 

The realisation of SBS in any device hinges on achieving sufficient Stokes amplification 



within the device length to be useful. Within the basic theory, where SBS is the only nonlinear 
process present, the amplification of the Stokes wave (i.e. the ratio of Stokes power at the output 
to injected Stokes power) is proportional to the power of the injected pump beam; the Stokes 
wave initially exhibits exponential growth until it starts to deplete the pump. In reality, however, 
neither this exponential growth, nor the predicted linear relationship between amplification and 
pump power, can be expected over the whole range of pump powers. The reason for this is 
nonlinear loss in the waveguide, where third-order processes - particularly two-photon absorp¬ 
tion (2PA) - impact the SBS-performance quite differently from fifth-order processes such as 
three-photon absorption (3PA) and 2PA-induced free carrier absorption (FCA). The latter is a 
major nonlinear loss mechanism in group-IV semiconductors such as silicon and germanium 
in important wavelength ranges. These materials are of interest because many promising ap¬ 
plications of SBS require the integration of SBS-gain in group-IV photonic circuits in order 
to benefit from CMOS-type fabrication techniques for mass production and integration with 
electronic circuitry. 

In this paper, we apply recent analytic work [18] to examine the limitations that nonlinear 
loss imposes on the gain for SBS. The related effect of nonlinear loss in Stimulated Raman 
Scattering has been investigated [19] although based on different approximations and with an 
emphasis on short-pulse propagation and pump depletion. In contrast, we present two important 
findings: firstly, that there is an upper bound for the output Stokes power that can be obtained 
with a specific waveguide design (i.e. for a specific geometry of the waveguide cross-section 
with a specific choice of materials) by amplifying a weak (e.g. thermal) initial Stokes signal. 
This bound, which involves both nonlinear loss coefficients and the SBS-gain parameter of the 
design, in particular represents an absolute upper bound for the output power of any continu¬ 
ous wave (CW) Brillouin-laser. Secondly, we find that the total amplification (ratio of output 
to input power) of an incident Stokes wave that can be obtained with a specific waveguide de¬ 
sign also has an upper bound. This maximally realisable amplification and the corresponding 
optimal operating conditions are determined by a d i m e n s i o n I e s s/z# « re of merit, which involves 
only the effective coefficients for the linear waveguide loss, the 2PA, the 2PA-induced FCA 
and the SBS-gain coefficient. Our figure of merit forms a valuable engineering tool for the de¬ 
sign of SBS-active waveguides since it provides a convenient way to compare the performance 
of different designs and yields good estimates for the waveguide length and pump power. We 
demonstrate this in the last section of the paper using two types of suspended silicon nanowires 
in backward SBS configuration as examples. 

2. Approach 

We choose to describe the problem within our coupled-mode framework [20] for SBS in an 
integrated waveguide, which is assumed to extend in the positive ^-direction. When focusing 
on strongly confining waveguides, as in this paper, it is appropriate to express the optical fields 
in terms of power carried by the individual waveguide modes, although other formulations are 
possible, e.g. based on field intensities as is common in bulk nonlinear optics. The equations 
of motion for the power P^ l \z) carried in the Stokes wave and the power P (1> (z) in the pump 


wave are 


sd z P w = (r-2j3-rP (2) )F (2) F (1) -aF (1) +Ci, 
<3 Z P (2) = -(15 + jP {2) )[P {2) } 2 - aF (2) + C 2 , 


(1) 

( 2 ) 


where the symbol s is used to describe both forward SBS (,v = +1) and backward SBS (s = — 1) 
within a common framework. Here T = f dr 2 u* • f is the overlap integral between the optical 
force distribution (involving both electrostriction and radiation pressure) as known from the 


literature [7,20]. The acoustic angular frequency Cl is defined by the acoustic dispersion relation 
and the optical wave number via the usual phase matching condition [7,20]. Later we will also 
provide theoretical estimates for the Brillouin line width Ail, which we compute from the 
overlap integral [20] between the acoustic eigenmodes and the dynamic viscosity of silicon and 
germanium [21]. 

The linear loss coefficient a, 2PA-coefficient p, and coefficient y for 2PA-induced FCA are 
overlap integrals of the optical eigenmodes and the respective dissipative nonlinear polarization 
currents as presented e.g. in Ref. [18]. 

The two large-signal correction terms 


Q = - ((3 + yP (1) +4yP (2) )[P (1) ] 2 , 

C 2 = - (2/3 +T + 4yP (2) + r pW)pWp( 2 ) 


(3) 

(4) 


can be neglected when assuming that the Stokes wave is weak compared to the pump wave, 
which we refer to as the small-signal approximation. Within this approximation, forward and 
backward SBS yield the same results for the power and amplification limits. We have assumed 
quasi-CW light and have absorbed the acoustic envelope into T via a local acoustic response 
approximation. This is justified if the optical modes vary on a length scale large compared to 
the acoustic decay length [20]; in particular, the optical loss must be sufficiently low [22]. 

Equations (1—4) are derived in detail in Ref. [18] and we restrict ourselves here to a brief 
description of the nonlinear loss terms. The terms involving p describe 2PA due to power 
carried in the respective mode itself (prefactor P) and due to the combination of both modes 
(prefactor 2/3). The nontrivial factor 2 arises because interference between the modes leads to 
a non-uniform intensity distribution; it appears commonly in nonlinear optics [19]. The terms 
involving y describe fifth-order loss (mainly 2PA-induced FCA) and as before appear as loss 
induced by either mode alone (prefactor y) and cross terms that describe the contribution of 
inter-mode interference to the 2PA-induced FCA (prefactor 4y). The nontrivial factor 4 arises 
in a similar fashion [18] as in the case of 2PA, but is less commonly encountered. 

For waveguides with eigenmodes that differ strongly from plane waves (e.g. silicon 
nanowires), the coefficients T,a,P and y depend sensitively on the geometry and eigenmode 
field distributions in the transversal plane. In the case of waveguides with plane wave-like eigen¬ 
modes (e.g. silica fibers), the optical powers can be replaced with local intensities and the pa¬ 
rameters with established parameters [1,4]. This is equivalent to normalising the eigenmode 
(i.e. plane wave) over a cross section of 1 m 2 within our framework. 

Finally, we note that inter-mode SBS can be described by the same equations at the expense 
of a more convoluted notation, leading to analogous results with modified parameters (see ap¬ 
pendix). 

3. Maximal output power 

We begin with the output power upper bound. From the large-signal Stokes equation 
[i.e. Eqs. (1,3)] we find 


sd z P w < [r-2j3-4yp( 1) -yP (2) ]P (1) P (2) . 


(5) 


In conjunction with the requirement sd z P * ')> 0 expressing Stokes amplification, this provides 
a necessary (though not tight) upper bound for the output power: 



F — 2/3 — 4y/4 2 ) ^ 

4 y < 4 


( 6 ) 





Fig. 1. Optimal pump power (left panel) and optimal waveguide length (right panel) as 
functions of the figure of merit [see Eq. (12)]; numerical results (black solid lines) in 
comparison to the approximate expressions Eq. (14) and Eq. (13) (red dashed lines) The 
optimal pump power is between £?/2 and 3? and the optimal waveguide length is of order 


where we have introduced the natural unit of power 

& = (r-2p)/y. (7) 

Any greater Stokes power level will necessarily decay with propagation and the maximum 
power limit can only be exceeded anywhere in the waveguide if it is exceeded at the input face. 
Notably, Eq. (6) is an absolute upper bound for the output power of any CW Brillouin-laser. 
This should be distinguished from the influence of linear loss, which mainly influences the 
laser threshold and the line width. Note also that the simple expression &/4 overestimates the 
maximal output power because of the term 4yp( 2) and further nonlinear loss terms, which were 
removed in arriving at the inequality Eq. (5). 

4. Maximal amplification and optimal design parameters 

The discussion of optimal design parameters and figures of merit is based on the small-signal 
equations, i.e. Eqs. (1,2) under the approximation C\ = C 2 = 0. As a result, we obtain upper 
bounds for both the pump wave and the Stokes wave power levels, and, therefore, for the mea¬ 
surable total Stokes amplification of a given waveguide. This can be concluded from the fact 
that the neglected terms Eqs. (3,4) are strictly negative. The small-signal equations can be 
solved analytically for the case j3 = 0 and starting from this solution, the case of j3 ^ 0 can be 
treated perturbatively [ 18]. In order to strip the problem of as many free parameters as possible, 
it is advisable to express all powers with respect to the natural power unit 3? introduced in 
Eq. (7) and to express all lengths with respect to the corresponding natural length unit 

A = 7 / (E — 2/3) 2 . 


(8) 


















Fig. 2. The maximally realisable Stokes amplification for FCA-dominated setups as a func¬ 
tion of the figure of merit & [see Eq. (12)] computed numerically (black solid line) and 
approximate fit according to Eq. (15) (red dashed line). Materials and waveguide designs 
with & < 1 cannot amplify an injected Stokes wave. 


We thus obtain an analytical expression that predicts the total Stokes amplification srf as a 
function of the waveguide length L and the injected pump power l 1 ' 2} (0) = Pq: 


srf(L,Po) =101og 10 
10 


pW(L) 

PM(0) 


dB 


tan 1 i/(l + t//)exp(2aL) — 1 — tan 1 1 // 


In 10 [_ y/aA - 

(1 + y/)exp(2aL) — 1 


— - In 
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¥ 


dB. 


(9) 


( 10 ) 


where t/r = a / (y/jj) is a measure of the relative strength of linear loss and 2PA-induced FCA. 
The factor 1/y/aA = (T — 2j3)/ v / ay in front of the first term of Eq. (10) is of particular in¬ 
terest, because it arises naturally from Eq. (1) by requiring that a positive P' 2i exists such that 
sd z pW > 0 (i.e. Stokes amplification at some positive pump power). Since all coefficients are 
positive, this is only possible if 


pi 2 ) = . 


2y 


(r — 2/J) ± y/(2/3 — T) 2 — 4ay 


> 0 ; 


r— 2/3 

2 v'oy 


>1. 


( 11 ) 


This is a necessary condition for any amplification of the Stokes wave. Therefore the quantity 


_ r-2/3 

” 2yay 


( 12 ) 


lends itself as a figure of merit for SBS in waveguide designs. 

Every jF corresponds to a different maximally realisable Stokes amplification ,fe/ !max - :i and, 
through (10), this amplification is only obtained for a specific (optimal) choice of waveguide 
length L' opli and pump power / J ' :opt - : '. We determined these values by numerically finding the 
maximum of ,e/ (L,Pq) for a wide range of & and show them as black solid lines in Fig. 1 and 
Fig. 2. Observe that the worse the figure of merit, the weaker the pump and shorter the wave¬ 
guide that delivers the optimal performance: since FCA is a fifth-order process while the SBS 






























Table 1. Table of nonlinear coefficients (T, j 6 , y), acoustic frequency £1/2k, Brillouin line 
width A£i/27t, SBS-figure of merit , maximally realisable Stokes amplification 
and the natural unit of power 3? (four times the max. SBS-laser output power in case of 
& > 1 ) for bulk silicon and germanium at wavelengths on the red and the blue side of 
the 2PA-threshold. We assumed a linear loss of a = 0.1 dB/cm and a carrier life time of 
10ns throughout. The quantities T and £2 were computed from literature expressions [20] 
assuming wave propagation along the [ 100 ]-direction and purely electrostrictive coupling 
based on literature photoelastic coefficients [24,25]. The Brillouin line width was com¬ 
puted as described in the main text using literature values for the dynamic viscosity [ 21 ]. 
The bulk 2PA-, 3PA- and FCA-coefficients are literature values [26-28] (annotated in the 
square brackets). As mentioned in the main text, this table lists bulk parameters, which are 
naturally expressed as power intensities (units of W/m 2 ) rather than powers; this is reflected 
in the units for si, 7 and 3? and p(°P l ), which differ from those in Table 2. Note that 
the photo-elastic tensor and the dynamic viscosity are available in the literature only for 
selected wavelengths, so we assumed them to be nondispersive, leading to a wavelength- 
independent prediction of the SBS-gain. 
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* three photon absorption is the leading multiphoton process at this wavelength 
t assuming 10 ns carrier life time 

t only a very rough lower bound since our fit [Eq. (15)] only covers the range & < 7 


gain and 2PA are third-order processes, trying to drive a poor system harder is counterproduc¬ 
tive. 

The numerically obtained relationships between ^ (in the range & > 1) and the performance 
and operation parameters can be represented by three phenomenological fit functions 


p(opt) f 

a (2 In J^) 0 ' 713 or 1 , 


(13) 

pi opt) f 

■('- 
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U> 

O 

y/&(&-!)+ 3.0 
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(15) 


These expressions are provided as a convenient way to predict the maximal amplification and 
optimal waveguide length and pump power in the process of designing an SBS-active wave¬ 
guide. Their good accuracy over a wide parameter range is demonstrated by Fig. 1 and Fig. 2, 
where they are shown as red dashed curves. 

As mentioned, the interpretation of Eq. (12) as a figure of merit is based on the fact that a 
material or waveguide design with JF < 1 cannot amplify a Stokes wave. For /# > 1, Stokes 






















Table 2. Table of SBS-resonance parameters, loss parameters, SBS-figure of merit, max¬ 
imally realisable Stokes-amplification and optimal operating conditions for two simple 
suspended nanowire designs (see Fig. 3) operated in backward SBS configuration in the 
2PA-regime. The effective waveguide coefficients were computed in analogy to and using 
material parameters listed in Table 1, yet using the specified optical and best-suited acoustic 
eigenmodes. The SBS-gain coefficient includes both electrostriction and radiation pressure 
terms as described in the literature [20], Furthermore, we show data for an experimentally 
studied geometry [9], where we adopted the SBS-parameters (gain, shift, linewidth), the 
linear loss and the carrier life time published for that particular structure and computed the 
nonlinear loss and from this &, and and extrapolated ,s/( max ). Optimal operation condi¬ 
tions cannot be provided, because this structure does not provide net Stokes amplification. 
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* geometry is scaled with the operating wavelength 
t computed assuming the carrier life T c 
t value taken from Ref. [9] 

o extrapolated using Eq. (15) slightly outside its intended parameter range 


amplification is possible, but the realisable amplification factor is limited to a number given 
by Eq. (15). For values 0 < & < 1, SBS generation of the Stokes wave occurs in pump-probe 
experiments with externally injected Stokes seeds, but with net attenuation along the wave¬ 
guide. However, the regime of & > 1 can still be reached by reducing the linear or fifth-order 
loss terms e.g. via free carrier extraction techniques. Finally, for & < 0, SBS is overcome by 
the competing two-photon absorption irrespective of pump power (since both SBS and 2PA 
are third-order processes) or waveguide length, and no SBS-allowed regime can be reached by 
adjusting a or y. 

5. Examples and Conclusions 

We demonstrate the utility of our figure of merit by discussing bulk silicon, bulk germanium, 
and three (two theoretical, one experimentally realised) silicon nanowire systems. In the case of 
bulk materials, the previously presented theory applies directly on replacing powers with power 
densities and effective waveguide loss and gain coefficients with the respective bulk material 
parameters. This is presented in Table 1 for the case of silicon and germanium at wavelengths 




















around the 2PA-threshold. Since this table contains numbers for SBS in bulk, only backward- 
SBS is possible and the SBS-gain is entirely defined by electrostriction and photoelasticity; 
radiation pressure does not contribute in Table 1. The SBS-gain and acoustic parameters were 
computed according to established methods for waveguides [7,20] using periodic boundary 
conditions in the transverse plane and literature values for the mechanical and photoelastic pa¬ 
rameters [23-25]. The acoustic line width was computed from the acoustic group velocity 
and the acoustic decay parameter « ph , which in turn was computed from the dynamic viscos¬ 
ity [20,21]. The nonlinear loss parameters are literature values [26-28]. It is clearly visible that 
bulk silicon and germanium are ill-suited for SBS on the blue side (shorter wavelengths, higher 
frequencies) of their respective 2PA-thresholds (« 2300nm and « 3650nm, respectively). This 
is mainly due to the combination of insufficient photoelastic coupling, high acoustic frequency 
(entailing high mechanical loss) and the strong fifth-order absorption term caused by 2PA- 
induced free carriers. However, we find that the fifth-order loss caused by 3PA allows for 
substantial figures of merit and should be negligible in practice. Therefore, SBS in bulk and 
quasi-bulk silicon and germanium seems infeasible in the near-IR and should be studied in the 
mid-IR. 

Furthermore, we present numerical calculations of figures of merit for two theoretical 
nanowire designs at two different wavelengths and one experimental silicon nanowire at 
1550nm. The two theoretical structures together with the relevant optical modes are depicted 
in Fig. 3; note that the second example is not operated with the fundamental optical mode. 
For both waveguides, we computed the backward-SBS coupling (involving both electrostric¬ 
tion and radiation pressure) with the lowest few acoustic modes and present the results for 
the lowest quasi-longitudinal mode. Again, we computed the SBS-parameters and effective 
nonlinear coefficients using established methods [7,18,20] and literature material parame¬ 
ters [21,23,24,26]. The results are shown in Table 2. For the theoretical structures, we have 
to choose a value for the linear loss, a quantity which is still improving with fabrication de¬ 
velopments. We present results for 1 dB/cm (as a typical current value) for calculations at both 
wavelengths (columns 2,3,5,6), and for 0.2 dB/cm as a state of the art value for the calculations 
at Ao =2000 nm only (columns 4,7). 

The results indicate that both designs can support SBS (i.e. & > 0) at both considered wave¬ 
lengths. However, neither of them is expected to provide net Stokes amplification on the blue 
side of the 2PA-threshold unless the linear loss is reduced considerably below 1 dB/cm. We 
find it noteworthy that our proposal for a waveguide based on the £ v -po lari zed eigenmode is 
predicted to provide higher SBS-gain than the previously studied £,’ v -polarized eigenmode. 

For the comparison with the experimental structure, we adopted the geometry, published 
forward SBS-parameters and linear loss of a nearly suspended silicon nanowire [9] and com¬ 
plemented this with our numerically computed nonlinear loss coefficients based on the carrier 
life time measured in that structure. The first observation is that we find a figure of merit just 
below 1 corresponding to an extrapolated total Stokes amplification of — 0.2dB in agreement 
with the published value of —0.1 dB. As the natural power unit we find P/ = 18mW, which is 
1.5dB less than the published 2PA-threshold of 25 mW; however, care should be taken when 
comparing the results, since the experimental 2PA-threshold gain was determined from Fig. 3a 
in Ref. [9], which depicts an on/off-gain, while we derived P/ in the context of optimal net 
Stokes amplification. 

In their paper, the authors of the experimental study mention that their device is close to net 
amplification. In fact, a reduction of the linear loss from the reported 2.6dB/cm to 1 dB/cm 
would lead to a & = 1.5, which would allow for a maximum net Stokes amplification of 2.7 dB 
realised for a waveguide of 37 mm and a pump power of 16mW. In contrast, a reduction of the 
carrier life time to 2 ns would lead to a similar ^=1.6, but the maximum Stokes amplification 
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Fig. 3. Waveguide geometries and effective indices of the optical modes that are studied 
for backward SBS in Table 2. The color-plot depicts the modulus of the modal electric 
field in arbitrary units, the black arrows indicate the in-plane electric field components. 
The waveguides consist of silicon in [100]-orientation; their cross-sections are shown as a 
thin black rectangle and are scaled with the light's vacuum wavelength A<). 


of 3.4dB would be achieved for a length of only 16 mm and a considerably higher pump power 
of 48 mW. The output powers of hypothetical CW SBS-laser using those two improved struc¬ 
tures would be limited to less than 9mW and 27 mW, respectively. This example illustrates how 
the results and parameter fits presented in this paper can help in designing the optimal length 
of samples and hint at the best route to improve performance. 

In summary, we have studied the impact of third-order and fifth-order loss on the process of 
SBS in semiconductor waveguides. First, we derived an upper bound to the output power of 
an SBS-based amplifier. Second, we introduced a figure of merit as a simple measure of the 
suitability of an SBS-design and provide simple fit functions that predict the SBS-performance 
of that structure as well as the optimal operating conditions. 
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covery Grant DP130100832, its Laureate Fellowship (Prof. Eggleton, FL120100029) program 
and the ARC Center of Excellence CUDOS (CE110001018). 


A. Appendix: Expressions for inter-mode coupling 

Within the main text, we restrict ourselves to intra-mode coupling with the goal of cleaner 
notation. In this appendix, we provide the corresponding results for the general case albeit at 
the expense of a great number of subscript indices. The governing equations (see Ref. [18]) for 
the general case of inter-mode coupling are: 
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Note that in contrast to intra-mode coupling, where the same loss coefficient (a, /3, 7 ) appears 
in several terms, these symmetries are now lifted because the coefficient in each loss term is 
associated with a different combination of eigenmodes. These are the power-related coefficients 
required in Eqs. (16,17) and have been derived from the amplitude-related coefficients derived 













in the Appendix of Ref. [18]: 
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Here, e 1 2| represents the eigenmode patterns, J 2 ^ 1 / 2 ) the corresponding power flux of the 
normalised eigenmode and E 2PA and E FCA are the nonlinear conductivities used to describe 
two-photon absorption (2PA) and 2PA-induced free carrier absorption (see Ref. [18]). 

The Stokes power limit follows from the full Eq. (16) exactly in the same way as described 
in the main text: 


p(i) < F 2 ^ 12 . (21) 
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The figure of merit follows from the small signal terms of Eq. (16) alone, again exactly 
as described in the main text for the case of intra-mode coupling [note that the /-coefficient 
appearing in the denominator is not necessarily the same as the one in Eq. (21)]: 
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